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Exact quantum master equation for a driven Brownian oscillator system is constructed via a 
Wigner phase-space Gaussian wave packet approach. The interplay between external field and 
dissipation leads to this system an effective field correction that enhances the polarization. This 
cooperative property is resulted from an effective bath response to the external field applied on the 
system. It is important in the low-frequency driving and intermediate bath memory region. We 
demonstrate this non-Markovian effect on the linear response and nonlinear dynamics and analyze 
the results together with analytical asymptotic expressions. 
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I. INTRODUCTION 

Brownian oscillator is the simplest and exactly solv- 
able system in quantum dissipation theory^! 2 ^^ Phys- 
ically it serves as an elementary model in diversi- 
fied fields, such as quantum optics, quantum statis- 
tics, spectroscopies, and chemical dynamics in condensed 
phases^LMiiaJiiSiiiii In this model, the harmonic 
oscillator (q) is coupled to a heat bath of harmonic os- 
cillators ({xj}) in a bilinear form. The total composite 
material Hamiltonian reads 2 ^^ 
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The effect of bath is completely characterized by the in- 
teraction bath spectral density function: 2 



mjCUj 



(2) 



Note that the Caldeira-Leggett model, Eq. (flaj) . is 
intrinsically perturbative, as the bath is assumed to 
be harmonic, i.e., h B = J2j[p"j/(^ m j) + m 3 U) 1 x ' l j V 2 ]> 
and linearly coupled to the system. However, it can 
be treated exactly without further assumption for an 
arbitrary system Hamiltonian, including anharmonic 
cases i£ i 15 i 16 i 17 One can obtain the Feynman- Vernon in- 
fluence functional path integral expression 2 ^ or its 
equivalent differential formalism in terms of hierarchi- 
cal equations of motio n 16 ' 17 for the reduced system dy- 
namics. Great simplification does exist for the Brownian 
oscillator where the reduced system of Eq. (|lb[) is also 
harmonic^ 19 ! 20 ! 21 ! 22 

The dynamics of Brownian oscillators can be studied 
in terms of the quantum Langevin equation. It starts 



with the Heisenberg equations of motion for all degrees 
of freedom, followed by using the formal solutions to the 
bath ones. The resulting quantum Langevin equation 
reads 2 

p{t) = -Am 2 q{t) - f dr 7 (t - t)p(t) + F{t - t Q ) 

-Mj(t-t )q(t ). (3) 



Here, j(t) = M 1 : {<q / rrijUJj) cos(wji) is the friction 

kernel function; F(t) = e ihBt / h (J2j Cj%j)e~ ihBt/H is the 
stochastic force operator, and its correlation function is 
related to the bath spectral density and the friction kernel 
functions via the well-established fluctuation-dissipation 
theoremjii 2 -^ Note that if the initial time is set to be 
to — > — oo, the last term in Eq. ([3]) is reduced to zero, due 
to the friction kernel function j(t — > oo) = 0. 

The Brownian oscillator system is more often stud- 
ied in terms of the quantum master equation (QME). 
The key quantity here is the reduced density oper- 
ator p(t) = ti B pr(t), i.e., the bath subspace trace 
of the total density operator. The exact QME to 
Brownian oscillators had been constructed via vari- 
ous approaches. These include the calculus-on-path- 
integral algorithm ; 18 i 19 the trace-over-bath on total- 
space Wigner-f unction method ; 20 ' 21 ' 22 the stochastic 
field 2 ^ and the stochastic Schrodinger wave function 2 ^ 
approaches. All these cited efforts had not considered 
the effect of time-dependent external fields. 

We have recently constructed an exact QME to a 
driven Brownian oscillator (DBO) system, with the as- 
sumption of initial thermal equilibrium before the time- 
dependent laser field action^ 2 ^ In this work, we will re- 
visit the theory with an arbitrary initial factorized state. 
A single-mode DBO system will be considered for clar- 
ity. It can be readily extended to multi-mode systems, 
together with the applications to some nonlinear spectro- 
scopic problems in condensed phases^ ' 26 ' 27 

The remainder of the paper is organized as follows. 
With the derivation detailed in Appendix^ we present 
the final results of QME-DBO theory in Sec.flll We 
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discuss its Markovian simplification in Appendix[B] and 
its asymptotic equilibrium property in Appendix[UJ The 
new feature of the theory is an additional effective field 
correction, resulting from the interplay between external 
field interaction and non-Markovian dissipation. We ana- 
lyze the effects of the field correction on the DBO system 
in Sec. lIIIi and finally conclude the paper in Sec. lIVI 



The frequency-domain resolution of Eq. ([8]) is then 

x(lo) = {Mp 2 -lu 2 -iujfa)]}- 1 . (10) 

Its real and imaginary parts are corresponding to the 
dispersion and absorption, respectively. The response 
function x(t) w iU be the key quantity in the following 
construction of the QME-DBO theory. 



II. THEORY 



B. Quantum master equation: Exact results 



Model and notations 



The total DBO Hamiltonian in the presence of an ex- 
ternal time-dependent field e(t) reads [cf. Eq. |T])] 

H T (t) = H M - qe(t) = H(t) + h B - c jXj . (4) 



Here, H(t) = Hq + H Icn — qe(t) is the sum of those 
in the system space, with the renormalization H rcn — 

q 2 Y,jCj/( 2m j u j) = \ M l{Q)l 2 - Thus > the system 
Hamiltonian in Eq. Q is 

H(t) = -^ + \m[Q 2 + 7 (0)]g 2 - qe(t). (5) 

The effects of bath interaction are described by the force- 
force correlation function: 



C B (t) = (F(t)F(0)) B = h[a(t) + ib(t)]. 



(6) 



Here, a(t) and b(t) denote the real and imaginary parts 
of C B (t)/h, respectively; (-) B = tr B (-pl q ) with p c B q = 
e~ /3/lB /tr B (e~ ,3 ' lB ) being the bath canonical equilibrium 
density operator, and (3 = l/(k&T) the inverse temper- 
ature. The imaginary part is the sine transformation of 
the spectral density (with a factor of —1/ir), and related 
to the friction kernel function as b(t) = 4f-7(i). The 
real part, a(t), depends on temperature and is related to 
the imaginary part, b(t), via the fluctuation-dissipation 
theorem.—^ 

For later use, we denote the response function of the 
Brownian oscillator system as 



x(t) = *<[g(t),g(0)]> M //L 



(7) 



Here, q(t) = e iHut/h qe -iH M t/h and ^.^ = tr M (-p^ q ) = 
tr M (-e~ / " fM )/tr M (e~ ,3ffM ), defined in the total system 
and bath composite material space. From its definition 
and Eq. ([3]) (setting to — there), we have 



X(t) + n 2 X (t) + / dr 7 (t - r)x(r) = 0, (8) 
Jo 

with x(0) = and x(0) = 1/M. Let 



m= dte^ X (t), Kw)= dte^j(t). (9) 
Jo Jo 



For convenience, we set to = hereafter, and assume 
the total system-bath composite state is initially of an 
uncorrelated form of pr(0) = p(0)p B q . Due to ergodic- 
ity the effects of initial states vanish completely in the 
course of time and the system relaxes to the thermal 
equilibrium. In Appendix[S] we derive the exact QME- 
DBO theory via the modified Yan-Mukamel method on 
the Gaussian wave packet (GWP) dynamics^ ' 25 ! 26 

The final results are summarized as follows. 
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-[H eS (t),p]~TZ(t)p, 



(11a) 



with 



H cS (t) = j M + \Mn 2 t q 2 - q[e(t) + 6e(t)}, (lib) 



and 



nt)p = ^h{p,p}]+^hhp}}-jhip,p}}- (He) 

The involving time-dependent functions are 



(12a) 
(12b) 



o2 _ X 2 (t) - X (*)*(*) 
* X 2 W - *(*)*(*) ' 

r = x{t) x {t) - mm 
1 x 2 (t)-mx(t) ' 

C t '= / drx(T)e(r;t), Ct = / drxW^r;*), (12c) 
Jo Jo 

Se(t) = f dT X e(t-T;t)e(r), (12d) 
Jo 



with 



£(t; t) = o(t) + / At' a(r - r')x e (r'; t), (13a) 



(13b) 



Xc (r; t) = M[x(t) + r t x(r) + n 2 X (r)], 



C. Comments 



The H c g(t) of Eq. (|llb[) . even in the absence of Se(t), 
is neither H(t) [Eq. ([5])] nor its Caldeira-Leggett coun- 
terpart, Ho — qe(t) [cf. Eq. (fib])] . The latter is recov- 
ered in the Markovian limit; see Appendix[B] The val- 
ues of x(0) = x(0) = 0, x(0) = 1/M, and X(0) = 
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FIG. 1: Real and imaginary parts of the linear response func- 
tions x(u>) [Eq. (|10[) ]. with the coupling bath strength V — 1 
and two values of spectral width (in unit of M = Q — H — 1): 
W — 2 (solid) and 5 (thin-solid). The dashed curves are the 
<5e-free approximation counterparts. 

-[ft 2 + 7 (0)]/M lead to fi£ =0 = n 2 + 7 (0) and r t=0 = 0; 
i.e., H e s(t) assumes H(t) [Eq. ([5|)] at t — before dissipa- 
tion occurs. When the solvent fluctuation is much slower 
than the system response dynamics, Eq. {8]| is reduced to 
x{t) + [ft 2 +7(0)]v(t) = 0; the relaxation is yet to occur. 
This implies that in the slow modulation limit the bath 
contributes only to the static reorganization energy. 

Temperature enters the QME through £ t ' and Ct via the 
real part of bath correlation function, a(t) = KeC B (t)/H; 
cf . Eq. (|12cp with Eq. (|13ap . The asymptotic behaviors of 
£j and Q are related with the thermal equilibrium vari- 
ances (T^ an d &qqi see Appendix[C] This guarantees the 
correct equilibrium property of the QME. 

The main feature in the present QME-DBO formal- 
ism is the inclusion of the effective local-field correction 
Se(t) [Eqs.CGd} and fl3bjl ]. It results from an inter- 
play between driving and non-Markovian dissipation. In 
the absence of the external driving field, the above QME 
formalism is the same as the previous theory; see, for 
example, Refs. flU andflfl 



III. DEMONSTRATIONS AND ANALYSIS 

A. The effective field correction 

The effective field correction Se(t) is purely of non- 
Markovian nature, (i) It vanishes in the Markovian 
white noise limit (cf. AppcndixlBj): (ii) For a general non- 
Markovian bath, if the field correction Se were neglected, 
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FIG. 2: The asymptotic Qt->oo and T t ^oo as functions of W 
for T =1 (solid), 0.75 (dash), 0.5 (dot), and 0.25 (dash-dot), 
respectively. Shown in the insets are the transient Q t and r t 
as functions of time, with the coupling bath strength V = 1 
and two values of spectral width, W — 2 (solid) and 5 (thin- 
solid). All are in unit of M — Q = h = 1. 



the QME would lead to the system response function 
x(cj) a Markovian form, [Mfflf^^ — uj 2 — iujTt^oc)}^ 1 ■ 
The exact result of Eq. (fit))) can only be recovered with 
the inclusion of Se(t) in the exact QME-DBO theory; 
(iii) Se(t) is most important in the low frequency driving 
and intermediate bath memory region; (iv) The effect of 
Se(t) is cooperative, leading to the enhancement of po- 
larization. 

Note that the asymptotic limit throughout this pa- 
per refers to the t > t c regime, where t c denotes the 
characteristic time of the system response function hav- 
ing x(t > tc) ~ 0. We adopt hereafter the internal 
unit of M = CI = % = k B = 1, and the Drude bath 
model, j(t > 0) = TWe~ wt , for the purpose of numer- 
ical demonstrations, where F and W denote the bath 
coupling strength and spectral width, respectively. 

Figure 1 compares the system response functions eval- 
uated with and without Se(t). The difference between 
them is considerable in the intermediate bath memory 
region, but vanishes in both W <C CI (narrow bandwidth) 
and W ^> il (Markovian) limits (thus not shown). 

Before further demonstrations, let us examine fit and 
r t , as these two quantities are involved in Se(t) [cf. 
Eq. (|13b[) ]. as well as in the Brownian oscillator part of 
the QME. Their initial values are O t=0 = yjti 2 + 7(0) 
and r t=0 = 0, as discussed in Sec. Ill CI Figure 2 depicts 
the asymptotic values of ttt^oc and Tt^^ varying with 
the bath spectral width W and coupling strength T, to- 
gether with the insets showing Q t and T t as functions of 
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time at T — 1, W = 2 and 5. The oscillations of fl t and 
r t vanish at large W. The oscillatory frequency is about 
that of the system response function. For the asymp- 
totic behavior, — + when W <C Q and — > T 
when W ^> £1, while Sloo — > O in both limits. 

Let us now turn to the local-field correction, 6e(t) of 
Eq. (|12dp . Note that <5e depends linearly on the external 
driving field e(t). We may therefore analyze the effect of 
the field correction with a continuous-wave (cw) drive, 
e(t) = e sin(Wft), as an arbitrary field is composed of cw 
fields. 

Note that the continuous-wave driven Brownian mo- 
tion had been widely studied on two-state/doublc-well 
systems , 28 ' 29 Rydberg atoms in a waveguide,- 3 ^ as well 
as dissipative harmonic oscillators , 31 ' 32 with both linear 
and quadratic dipole moments. Although the additive 
term of an effect driving (that would lead to a form of 
Se) was mentioned, see for example, the Appendix B of 
Ref. [321 . its influence was not studied. And with the fur- 
ther Markovian approximation made there, the effective 
field correction vanishes and goes unnoticed, as consis- 
tent with the limiting case of the present work; see Fig. [3] 
below and the discussion follows. 

The asymptotic behavior of 8e(t) associating with the 
cw-field can be obtained as [cf. Eqs. (|12dj) and (|13bl) and 
noting that X*(uj) = x(— u))] 

Se(t) -» e Im[e^*x*(w/)] = re sm(u f t - 9), (14a) 

with 

/>oo 

Xe(uf) = / dTe^f T Xe(r;t - 00) = re ie . (14b) 
Jo 

In other words, x e (a>/) defines the relative magnitude r 
and the retarded phase 9 of the asymptotic 8e(t) with re- 
spective to the original cw driving field. Figure 3 demon- 
strates the r and 9, for their dependence on the bath 
spectral width W, the coupling strength T, and the cw- 
field frequency u) / . 

The relative magnitude r of Se t -*oo > as shown in the up- 
per panel of Fig. [31 exhibits turnover behaviors, as func- 
tion of bath spectral width W. On the other hand, it de- 
creases monotonically as the external cw-field frequency 
ujf increases (see the inset of upper panel), due to the 
cancellation for the highly oscillatory e(t). As the func- 
tion r(W) is concerned, both its peak and width increase 
rather rapidly as the system-bath coupling strength T 
increases. We may refer the case being of experimental 
significance if r = \5e/eo\ > 0.2. With the moderately 
strong value of Y = 1, the field correction would take a 
significant effect within 0.6f2 < W < 7f2, which often 
covers the experimentally important region. 

The retarded phase 9 of <5et^oo, as shown in the 
lower panel of Fig. [31 decreases from ir/2 at long mem- 
ory (W <C f2) regime, to in the Markovian (W 3> fi) 
limit; meanwhile it increases from to 7r/2 as the driv- 
ing field frequency increases; see the inset in the lower 
panel. These observations are consistent with the pic- 
ture that the retarded phase reflects the relative bath 
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FIG. 3: The amplitude scaling factor r (upper) and re- 
tarded phase 9 (lower), which characterize the asymptotic 
6e(i) —* reosm(u>ft — 9), as functions of W, with u)f = 1 
and T =1 (solid), 0.75 (dash), 0.5 (dot), and 0.25 (dash-dot), 
respectively. Shown in the insets are r and 9 as functions of 
the cw driving field frequency uif, evaluated with T = 1 and 
W = 2. All are in unit of M = Q, = h = 1. 



response time with respect to the inverse frequency of 
driving. 

We have thus clearly demonstrated that the effective 
field correction Se(t) results from the interplay between 
the driving and non-Markovian dissipation. Remarkably, 
the observed retarded phases are all of 9 £ [0, n/2] . Thus, 
the effective field correction is of cooperative effect that 
enhances the influence of driving. The physical origin of 
this enhancement is due to the bath-induced disperse of 
5q that leads to an effective larger polarization. Neglect 
of Se will underestimate the influence of driving such as 
the field-induced polarizations to be demonstrated be- 
low. 



B. Driven Brownian oscillator dynamics 

Consider now the dynamics of system in the presence 
of an external field. Let us start with the asymptotic 
behaviors under the cw-field drive. The mean values of 
the coordinate and momentum can be obtained as [cf. 
Eq. (S)] 

q t —►eoMe*'''*** ("/)], ( 15a ) 
Pt — » e Muj f Bje[e iUft X*^f)}- (15b) 
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FIG. 4: The evolution of qt (upper) and pt (lower) under the 
cw-field e(t) = eosin(o;/t), with eo = 1 and u)f = 1 (solid), 
0.5 (dash), and 2 (dot), respectively. The bath parameters 
are V = 1 and W = 2. All are in unit of M = H = H = 1. 
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FIG. 5: The evolution of the phase-space variances, <j qq (t), 
<J P q(t), and a pp (t), at temperature T = 1 (solid), 0.5 (dash), 
and 0.1 (dot). The bath parameters are V = 1 and W — 2. 
All are in unit of M = fi = 7i = ft B = 1. 



The system energy, in terms of (Hq) = ^jj(p 2 ) + 
^Mfl 2 (q 2 ), has the asymptotic expression of 

+ \M(u} n 2 )e 2 {Me iu,t r(uf)}} 2 - (16) 

The three terms in the right-hand-side of the above 
equation describe, respectively, the thermal average 
value, the cw-field induced static polarization with the 
system elongation of 5q = \eox(u>f)\, and the residue 
oscillatory breathing. In the absence of dissipation, 
x(u>) = [M(n 2 — cj 2 )] -1 , both the static polarization 
and the oscillatory breathing terms diverge at u)f = Q. 
This is consistent with the classical energy climbing pic- 
ture of a cw-field driven oscillator. In the presence of 
dissipation, the oscillatory breathing is asymptotically 
quenched at resonance, tOf = Q, while the static po- 
larization contribution remains finite. Both the field- 
induced energy climbing and breathing into and out of 
the system are compensated by dissipation. It is also no- 
ticed that the resonance condition is closely related to 
the choice of reduced system Hamiltonian. For exam- 
ple, if the reduced system energy is examined in terms of 
(H(t = 0)) = 2li(p 2 ) + ^Mfl 2 =0 (q 2 ), similar expression 
of Eq. (fTTJl) holds, but with the resonance condition being 
replaced by uif = f2t=o = \J& 2 + 7(0). 

Demonstrated below will be the time evolution of DBO 
systems, with both a ciy-field (cf. Figs.0HS]) and a half- 



cycle of the field (cf. Fig. [7]). The results will be re- 
ported in terms of the first- and second-order cumu- 
lants, {q t ,Pt} and {a qq (t), a pq (t), cr pp (t)}, that fully char- 
acterize the GWP of the DBO system (cf. AppendixlA]) . 
and also in terms of (Hq) or polarization that is exper- 
imentally observable. Note that {qt,Pt} do not depend 
on temperature, while {a qq (t), a pq (t), a pp (t)} do not de- 
pend on the driving fields; and their equations of motion, 
Eq. (|A6|) and Eq. (|A8|) . are independent. 

For the DBO system in the presence of a cw-field 
e(t > 0) = eo sin(u)ft), we set the initial GWP to have the 
first-order cumulants of {qt=Q,pt=a\ = {—2,0} and the 
second-order cumulants the same as their thermal equi- 
librium values. This mimics the scenario of ultra-fast 
UV/visible laser pulse excitation of an electronic ground 
state Gaussian wave packet onto a linearly displaced ex- 
cited surface, immediately followed by the cw infrared 
field. The wave packet dynamics on the excited surface 
under the influence of the ciu-field is then studied. 

Figures 4 depicts the evolution of q t and pt- Their 
asymptotic behaviors, as described by Eq. (TT5l) . oscillate 
with the cw-field frequency u)t. 

Figures 5 demonstrates the evolution of cr qq , a pq , and 
a pp . Although they were initially of the equilibrium val- 
ues and they are independent of the external field after- 
ward, the second-order cumulants do oscillate and take 
time back to equilibrium. This is because there was no 
excited state wave packet initially and the dissipation 
considered here only starts when there is a real occu- 
pation upon excitation. Therefore the variances of the 



FIG. 6: The evolution of the mean energy (Ho), driven by the 
cu>-neld with eo = 1 and uif =0.5 (shifted by +2), 1 (shifted 
by +1), and 2 (no shift). The bath parameters are V = 1 
and W = 2 at temperature T = 1. The dashed curves are 
the <5e-free approximation counterparts. All are in unit of 

m = n = n = k B = i. 



1 2 

t/(27t) 

FIG. 7: The evolution of (8H ) = (Ho) - (-ffo)eq, driven by a 
half cycle of the cm-Held with eo = 1 and uif =0.5 on an initial 
thermal equilibrium state at temperature T = 1. The bath 
parameters are W = 2, and F = 1 (solid) and 0.5 (thin-solid), 
respectively. The dashed curves are the <5e-free approximation 
counterparts. All are in unit of M = O = h = fca = 1. 



pumped GWP, as depicted in Fig.[5j do experience dy- 
namics solvation. 

Figure 6 shows the evolution of the system energy, 
(H Q ) = 2Xr(p 2 ) + ±MQ?(q 2 ). The asymptotic behav- 
ior had been described in Eq. (fl"(5)) and analyzed there- 
after. In the non-resonance condition, u>f ^ f2, the en- 
ergy flows in and out with the cw-field frequency Uf. 
This energy breathing is quenched at resonance, uif = Q. 
Shown in the dashed curves are the 5e-iree counterparts 
for comparison. Their difference is proportional to the 
field strength, and also depends on the difference between 
the exact x(uif) an d its approximation; cf. Eq. (|I6[) and 
Flg.ffl 

Consider now the DBO dynamics driven by a half- 
cycle of the sine-function field. The system is assumed 
to be initially at the thermal equilibrium. Figure [7] re- 
ports the resulting dynamics in terms of the polariza- 
tion/magnetization related quantity, (5Hq) = (Hq) — 
(iZo)eq. Again it shows that the effect of 5e is coop- 
erative, and the neglect of this correlated driving and 
non-Markovian dissipation contribution will significantly 
underestimate the field-induced polarization. 



and non-Markovian dissipation only appears as a field 
correction on the coherent motion. In a general anhar- 
monic system, this correlated effect is expected to have 
much rich implications on the dissipative dynamics even 
in the weak system-bath coupling regime, see for exam- 
ple, Refs. HU and [34]- The exact non-Markovian dissipa- 
tive dynamics of general systems under arbitrary external 
fields can be studied via the hierarchical equations of mo- 
tion approach i 17 ' 35 i 36 i 37 i 38 i 39 These include the quantum 
transport systems, such as the transient current under 
the time-dependent bias voltage reported recently ] 40 ! 41 
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IV. SUMMARY 

This work presents the exact QME-DBO theory, for 
which the time-dependent external field interaction is 
considered explicitly. It shows explicitly that the cor- 
related driving and dissipation result in an effective field 
correction that enhances the polarization. This coopera- 
tive property is intrinsically non-Markovian and reflects 
the effective bath response to external field applied to 
the system. It vanishes in both the Markovian and slow 
modulation bath limits, but it is important in the inter- 
mediate bath memory and low-frequency driving region. 

For the DBO system, the correlated effect of driving 



APPENDIX A: THE WIGNER PHASE-SPACE 
GAUSSIAN WAVE PACKET APPROACH 

The QME-DBO theory [Eqs. CP to (I3J)] can be de- 
rived via the modified Yan-Mukamel method on the 
Gaussian wave packet (GWP) dynamics as follows! 3 ' 25 ! 26 
The solution to the Langevin equation of the DBO sys- 
tem reads^ 5 . 



9(0) 
p(0) 



T(t) 

+ / drT(i-r) 
Jo 





e(r) + F(r) 



, (Al) 
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with T(t) a 2 x 2 matrix: 
T(t) 



M 2 x (i) Mx(t) 



(A2) 



The Langevin equation in the time-local prescription is 
then 



Pit) 



A(t) 

Here, e c g(t) = e(t) + Se(t) and 





e e s(t) + F eS (t) 



F eS (t)=F(t) + / dT X e(t-T;t)F{T), 



1/M 
-A/^ -T t 



A(t) = T^T" 1 ^) = 



(A3) 

(A4) 
(A5) 



with fl 2 , r t , and 5e(t) being given in Eqs. (fl2)) . 

We can then obtain the equations of motion for all the 
cumulants of GWP^ 



Qi 

Pt 



A(t) 



Qi 

Pi 





e e ff(i) 



W t = A(t)W t + W t A T (t) 

(F eS (t)Sq(t)) 



Re 



{F eB (t)5q(t)) 2{F cfi (t)5p(t)) 



(A6) 



(A7) 



Here, W t = [ ] denotes the phase-space vari- 

ances of GWP, <5g(£) = <?(i) — <Zt, 8p{t) = p(t) — p t , and 
( • } = tr M [ • pt(0)]. We also use the initial factorized 
state p T (0) = piO)p B q - From Eqs. flll]) and CM}, we 
have 



F eff (t) 



W) 

6p(t) 



= / dr 

o Jo 

t 



X(t - t) 
Mx(t - r) 



{C B (t-T) 



+ ^dr'x £ (i-r';t)C B (r'-r)}. 



Here, C B (t) is the bath correlation function of Eq. ©. 
Comparing with Eqs. (|12cp and (|13a[) . it is easy to find 
that Re(F eS (t)5q(t)) = Kt and Re(F cS {t)Sp(t)) = 
MK't- Substituting those into Eq. (|A7|) . we obtain 



W t = A(t)W t + W t A T {t) + h 



Ct 



(A8) 



Now using the standard Yan-Mukamel method^ we 
arrive at the QME-DBO in the Wigner representation: 



d t pw 



y q ii +Am ^-^Yp 



pw 



d 



,d 2 



T t —p + MhC t — + h(. 



d 2 



dp^ dp 2 
This is equivalent to Eq. (TlTj) . 



dpdq 



Pw- (A9) 



APPENDIX B: THE MARKOVIAN LIMIT 

In the Markovian white noise limit, j(t) = 2TS(t) or 
r y{uj) — r, the system response function [Eq. (110|) ] reduces 
to 



X(* > 0) = (MO)" 1 sin(m) e - rt/2 , 
where Cl = (Q 2 - T 2 /4) 1 / 2 . It results in 

n t = n, r t = r, Xe (r;t)=o. 



(Bl) 



(B2) 



Thus the effective field correction [Eq. (|12d[) ] vanishes 
in this limit. It is well known that cr^ diverges in 
the Markovian white noise limit without further the 
high-temperature assumption. 1,2 Same situation applies 
also to Ct and Ct [ci. Eq. (|C1[) ]. We have therefore the 
high-temperature limit of ait) = Mk-B,T^(t)/K, so that 
Ct' = Tk B T/h and Ct = 0. Thus, Eq. (JTTJ is reduced to 



Pit) = ~[H Q -qe(t),p(t)] - l ^[q, {p,p(t)}] 



MTk B T 



fa fa pit)}}- 



(B3) 



This is the Caldeira-Leggett equation^ which has also 
been extended to anharmonic systems with an arbitrary 
system-bath coupling mode Q a o 42 ' 43 



Pit) 



[H - £e(t),p(t)] 

Q,[Q,pit)}}- 



h 

Tk B T 



r 

2h? 



[Q,{[H ,Q},p(t)}} 



K 2 



(B4) 



Note that the system Hamiltonian involved here is H in 
the Caldeira-Leggett form. 



APPENDIX C: EQUILIBRIUM PROPERTY 

In this appendix we will show that Ct and C* have the 
following asymptotic expressions, 

H&^Ttog/M, hC t ^MS%o%-o%/M. (CI) 

Here, t c is the characteristic time of the Brownian re- 
sponse function such that xit > tc) ~ 0. Thus, Eq. (|A8|) 
behaves asymptotically as 



W 



t>tc 



A(t)(W t - W cq ) + (W t - W cq )A T (t). (C2) 



This guarantees the correct equilibrium property of 
QME. The involving equilibrium phase-space variances 



are = and^i 3 - 



qq 



PP 



Im[x(w)] 



TV 

M 2 h 

TT 



e -l3hu ' 
(jj 2 Im[x(o;)] 



(C3a) 
(C3b) 
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To prove Eg. (|C1[ ), let us recast Eg. (|12c[) with Eqs. (jT3|) Note that a(t) = Ke[C B (t)]/h is an even function and 
as follows. related to the bath spectral density via the fluctuation- 



■M\x(t)+T tX (t)] f dTa(t-T) X (r) 
Jo 

+ M x (t) { dra(t-r)x(r) 



dissipation theorem)^ 



«(*) = i J J. (c ** + e -*-t) . (C5 ) 



-M / dr / dr'x(r)a(T-r')x(r') 
Jo Jo 

+ Am 2 t I dr / dr' X (r)a(r - t')x(t'), (C4a) 
Jo Jo 

C t ' = MQjx(t) / dro(t-T)x(r) Consider now the t > t c regime wh ere x ft > t c ) « 

Jo 0, the first two terms of both Eqs. (|C4a|) and (]C4b|) 

/"* can be neglected, and every upper limit of the inte- 

+ Mx(t) j &ra(t-T)x{T) grations can be set to oo. Then substituting Eg. (fC5|) 

t t into Eq. (IC4[) . performing Fourier transforms, and us- 

+ MT t dr At 1 X {r)a{T -t')x{t'). (C4b) ing the identity Im[xM] = J(cj)|xH| 2 with noticing 

Jo Jo J{lu) — MwRe[7(w)], we finally obtain Eq. JUT}. 
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